Theoretical model of structure-dependent conductance crossover in disordered carbon 
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We analyze the effects of sp^ / sp^ bond-aspect ratio on the transport properties of amorphous 
carbon quasi- ID structures where structural disorder varies in a very non-linear manner with the 
effective bandgap. Using a tight-binding approach the calculated electron transmission showed a 
high probability over a wide region around the Fermi-level for sp'^-rich carbon and also distinct 
peaks close to the band edges for sp''-rich carbon structures. This model shows a sharp rise of the 
structure resistance with the increase of sp'^C % followed by saturation in the wide bandgap regime 
for carbon superlattice-like structures and suggests the tuneable characteristic time of carbon-based 
devices. 
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Introduction: Amorphous carbon (a — C) thin films 
have long been considered as very important materials 
for understanding the fundamentals of electronic struc- 
ture :lh(i mechanical properties 043- Depending 
on the preparation method, a very diverse nature of mi- 
crostructure can be formed ranging from diamond-like to 
graphite-like carbon. Significant efforts have been made 
to explain the optical properties of a — C films based 
on the estimated band gap and degree of disorder [Tol - 
[isj . There are also numerous interesting results showing 
a structure-dependent conductivity crossover in a — C 
films [T9l - [2^ . However, a rigorous theoretical model of 
electrical transport of a — C films related to the carbon 
microstructure is yet to be developed that can explain 
the experimentally observed highly non-linear variation 
of the conductivity with sp^ bond concentration [l^, [l^ . 
In carbon the electronic properties can be strongly in- 
fluenced by bond disorder which depends on the sp'^/sp^ 
bond aspect ratio and on the sp^ cluster size. Although 
there is no well defined structure for a — C films they 
have been described in many studies as a mixture of 



sp and sp"^ hybridized bonds [lOHlSl Il8j . A number 
of spectroscopic techniques, for example Raman, opti- 
cal absorption and photoelectron spectroscopy, have con- 
firmed this assertion [2j, [13, [31 ■ It was found that the 
size of the sp'^ — C clusters greatly influenced the band- 
gap [J- The optical absorption edge was described by a 
Gaussian plot based on the normal distribution of clus- 
ter size [23] . Films with low sp^ concentrations have also 
been described as quasi-lD polymer chains (e.g. trans- 
polyacetylene) d [iMl [H, [Ij- The effect of ID fil- 
amentary channels on electronic transport has been ob- 
served in low-dimensional a — C films and related devices 
[25j . However, no significant theoretical studies have yet 
been undertaken to investigate the tunneling properties 
unlike in 1_D molecular structures QMi- Based on previous 
studies we believe that a — C films can be described ef- 
fectively as a distribution of periodical alternation of the 
hopping energy and constant energy term. Starting from 
previous experimental claims for the variation of disorder 



with the energy gap [T^, [iJ, ) we establish the trend 
of localization length in a wide range of carbon films 
based on the calculated transmission coefficients T{E) 
and local density of states (LDOS). In this Letter we 
propose a Gaussian disorder analogous to previous stud- 
ies [ll[il[i3 but in a different manner which is directly 
connected to the microstructure of a — C films. Since 
the nature of topological disorder has not been clearly 
understood from previous works we briefly discuss this 
effect on electron transmission at the end of this Letter. 

Proposed Structure: We develop a quasi ID- 
dimensional a — C superstructure, which can be rep- 
resented by a network of narrow nano-ribbons with a 
length of 87 A for the zigzag and 75 A for the arm- 
chair directions (Fig. 1(a)). This structure is a quasi- 
one-dimensional representation of disordered a — C, for 
similar examples see 0, [27,, J8] which include some of the 
essential features specifically sp^/sp^ clusters and bond 
angle distortions. There are 4 segments of the sp^ struc- 
ture with widths ranging from 1 to 10 sites depending on 
the phase percentage. The superstructure is a mix of sat- 
urated and conjugated {sp^ — sp'^) areas, corresponding 
to a and tt bonds. The percentage of sp^ carbon struc- 
tures determining the average sp^ cluster size [2] can be 
promoted, for example by nitrogen doping. The effects 
of the lateral dimensions can be taken into account by 
considering quantization of the wave vector in this direc- 
tion and modifying the hopping terms in the longitudinal 
direction. 

Methodology: We use a tight-binding Hamiltonian de- 
scribing electrons confined in g — ID a — C ribbons as 



H 



), (1) 



where e„ is an on-site energy, which can represent an 
atomic energy as well as an external potential, cjj and 
c„ are operators of creation and annihilation of electrons. 
t„-i_i_„ is a hopping term between site n and n+l, which 
takes into account nearest-neighbor hopping that can be 
different from site to site reflecting the structural change. 
We seek the solution of the Schrodinger equation 
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Elpn = enlpn - tn+l.ni'n+1 " tn-l,n'4'n-l (2) 



on the tight-binding lattice with the known eigenvalue 
E = ei ~ 2t cos ki , where the unit cell length is fixed to 
1. ti and 6/ are the hopping term and on-site energy, 
respectively, corresponding to the left lead, which can 
be interpreted as the Fermi energy (Ef), and fc; is an in- 
coming wavevector. The boundary conditions imply that 
there is only an outgoing plane wave at the right lead, 
while at the other sites in the tight-binging lattice there 
are incident and reflected plane waves. Starting from a 
unique outgoing wave on the right lead ipn = e'*^''", we 
backpropagate this solution to obtain ipi and ?Ao using 
Eq. We get the incoming wave amplitude A and 

define T{E) = f^ffl^, (outgoing wavevector, kr). We 
calculate the transmission coefficient as a function of in- 
cident electron energy, where the transmission coefficient 
is the ratio of the outgoing to incoming probability cur- 
rent (since k vectors may be different for incoming and 
outgoing waves). Thus, the problem is different from that 
of finding the eigenvalues as the transmission coefficient 
is evaluated as a function of energy. The resonance en- 
ergies determined from transmission maxima correspond 
to eigen-energies of the system. 

The Carbon model: In this model a Hamiltonian with 
alternating values of the resonance integral produces con- 
duction and valence bands as well as a band gap [1^ . For 
the saturated region we introduce the resonance integral 
alternation by the hopping terms ij^^j^ and f^ax- the 
presence of the conjugated compound the effective gap 
decreases because of the electron resonance at energies 
lying inside the sp'^ phase band gap. The hopping term 

corresponding to the sp^ phase experiences distortion 
due to the structural disorder, which reflects the differ- 
ence in the bond length via the deformation potential. 
Topological disorder corresponding to the difference in 
the cluster size affects the resonance conditions via the 
boundary conditions for quasi-bound states [l^ . The dis- 
order of the sp^ phase is neglected due to its minimal ef- 
fect on the resonant electron transmission within the gap 
0. The broadening effects are mostly due to disorder in 
the well. As a matter of fact the resonant states have a 
larger amplitude inside the wells in this range of energies. 
The leads on the right and left sides are taken as ordered 
narrow ID graphite-type strips. The incoming electron 
energies are given with respect to the tt electron (sp^ — C) 
on-site energy 6-^=0 eV. Therefore, the position of the 
valence and conduction sp^ phase is not symmetric due to 
the different on-site energies for tt and a electrons and has 
a shift corresponding to the difference in these energies, 
as shown in Fig. 1(b) The hopping between sp^ and sp^ 



hybridized C-atoms is extremely significant due to the 
non-planar geometry of the sp^ structure [7] . We use the 
following tight-binding parameters: sp^ on-site energy 
Ecr = —0.9 eV, sp^—sp^ transfer integral: ^2-2 = 3.40 eV, 
sp^ — sp^ transfer integral tc—„ = 1.37 eV , sp^ — sp^ 
transfer integr als = 1-10 and C""" = 4.30 eV 
[26I [29j . Fig. 1(b) shows the levels of hopping terms 
along with on-site energies for atoms in sp^ (cr) and sp^ 
(tt) structures. 
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FIG. 1. (a) Microstructure of a — C consists of a and tt bonded 
carbon atoms, (b) The band diagram shows sp^ hopping term 
alternation (red and blue lines online) and sp^ hopping term 
level. The longer dashed line corresponds to sp^ on-site en- 
ergy, which is taken as the zero level, the shorter dashed line 
corresponds to sp^ on-site energy, (c) T{E) vs. incident elec- 
tron energy (E) showed a number of quasi-bound states of 
mainly conjugated (78 % sp^) (dashed curve, red on-line) 
compared to pure saturated {sp'^) structure (solid curve). 



Fig. 



1(c) shows T{E) of a pure saturated structure 
and an ordered structure with 78 % of sp-^ bonding con- 
centration where the zero level corresponds to sp^ on-site 
energy. The numerically calculated band gap and band 
widths of the saturated structure are E^ — 6.4 eV and 
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2.2 eV , respectively. For the mixed structure. 



there are several narrow peaks within the sp'^ band gap 
corresponding to quasi-bound states of the sp^ "wells". 
The peaks are very sharp with a high peak-to-valley ra- 
tio. In the case of a symmetric structure, T{E) always 
reaches 1 [1^. When the concentration of sp^ structure 
is high, the peaks are close to the Ef. On the other 
side of the concentration limit (low sp^ phase) , the peaks 
are close to the band edges of sp^ structure. Fig. [UJa) 
shows the first resonance peak position E^es (from the 
zero energy level) as a function of sp^ phase percentage, 
which has an inverse square dependence on the percent- 
age. Such a dependence reflects a linear increase of the 
average sp^ cluster size and a quadratic decrease of the 
energy associated with the resonance. 

Effect of Disorder: There is a significant effect of 
structural sp^ disorder on the transport properties of 
the sp'^ — sp^ carbon system. Disorder is an inherent 
micro-structural property of a — C films. It has been sug- 
gested that the density of the bonding and anti-bonding 
TT states could be represented with a Gaussian distribu- 
tion [1^ [2^. In particular the Urbach energy was 
described as a function of the width of the tt DOS peak 
at Ef although the link between disorder and the broad- 
ening parameter was not clear. The disorder parameter 
6, which takes account of bond length (angle) distortions, 
was calculated previously [Tsj . Micro-structural disorder 
is broadly separated into structural and topological dis- 
order although the later only applies for a large sp^ clus- 
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ters. In this Letter a typical hopping disorder parameter 
S for a Gaussian distribution of the hopping term is 
shown in Fig. [2l[b) for 3 cases, A, B and C as a func- 
tion of Et^, which covers the major possible behaviors of 
disorder in carbon structures. Models B and C are sug- 
gestions for the behavior of the disorder parameter based 
on the assumption that the disorder parameter can either 
increase or decrease with the corresponding to either 
sp^ rich or sp"^ rich structures. The values of the dis- 
order parameter for model A S can be extracted from a 
number of experiments (e.g., Raman G-peak linewidth) 
[3 113, [13 and we believe it is the result of two competing 
processes (B and C). On this basis, we have determined 
the behavior of the disorder parameter as a function of 
We therefore focus on model A as it corresponds 
to the physical nature of a — C systems. A number of 
sp^ — C percentages were chosen to cover a wide range. 
The corresponding Ej^ values (and hence corresponding 
values for S) were determined based on the relationship 
between sp^ — C % and Et^. We found that the values 
of E'tt were very close to the values of Eres- Here 2Et^ 
is the peak-to-peak energy spacing of bonding and anti- 
bonding DOS above and below the Ef [l^. For case A, 
the disorder reaches its maximum at about E^^ = 1.9 eV, 
followed by a slight decrease. This is an attribute of the 
structural relaxation, which occurs (also for case B) when 
sp^ clusters become relatively large [2, 3]. For case C the 
disorder increases sharply with i?^ corresponding to the 
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FIG. 2. (a) First resonance peak from the zero level -Eres 
as a function of sp^ phase percentage, (b) Hopping disorder 
parameter 5 of the sp^ region for cases A, B and C as a func- 
tion of E-n , where 2_E^ is the peak-to-peak energy spacing of 
bonding and antibonding states. 



rise of sp^ phase content. In our further calculations we 
associate Et^ with the energy of the first (from zero) res- 
onance peak, Eres ■ The proposed non-linear graph (case 
A) showed a significantly different trend of T{E) and 
localization length compared with the two other cases 
where a continuous decrease (B) or increase (C) of 6 with 
Ejr is considered. 

In Fig. ^ (a) we show T{E) for 73 % sp^ bond concen- 
tration, which corresponds to 2Eres = 1-1 eV and Gaus- 
sian disorder ^ = 0.10 eV^ averaging the values over a set 
of 500 runs. Other cases are shown in Fig. [3] (b) 42 %, 
2Er.es = 1.7 eV, S = 0.19 eV, (c) 24 %, 2Eres = 2.2 eV, 
S = 0.26 eV and (d) 7.8 %, 2Eres = 3.7 eV, S = 0.29 eV. 
With the decrease of sp^ bond concentration and increase 
of disorder in model A, we observe an increase in the 
2Eres value (the effective band gap), broadening of the 
peak width and a decrease in peak-to-valley ratio. With 
further decrease of sp^ bond concentration, the resonance 
peaks move closer to the sp'^ band gap and finally disap- 
pear. Models A and C show significant broadening of the 
resonant peaks compared to case B for low sp^ % since 
disorder is large (see Fig. ^c) and (d), insets) although 
it is very similar for all cases for high sp^ %. At large 
sp'^ percentages the transmission through the structure 
is high therefore structural effects only slightly influence 
the transmission coefficient. When the sp^ percentage is 
low, the transmission is also low hence even slight struc- 
tural changes have a large effect on T{E). 

The disorder induces localization of the wave function 
within the structure and the localization length is ex- 
pressed as Lioc{Eres) = ~\^p^: — ) I where L is the length 
between leads used in the Landauer formula for calculat- 
ing the resistance [13]. At small energies Li^c is longer 
for case A and C than for case B because the structure 
is more transparent at resonant energies for lower disor- 
der (Figllja), lower curves). The localization length de- 
creases sharply for all 3 cases till Eres reaches the value 
of about 1 eV with a further saturation showing an ex- 
ponential behavior. The resistance behaves oppositely to 
the localization length (Fig|4l^a), upper curves). At small 
energies the resistance is smallest for case C and largest 
for case B following the distortion parameter dependence. 
The structural resistance increases till Eres = 2 eV fol- 
lowed by a slight decrease at higher energies in the natu- 
ral logarithm scale. This value is related to the mobility 
edge- as the disorder parameter increases, the transmis- 
sion maxima shift away from Ef towards the mobility 
edge. The lifetime of resonant states therefore decreases 
(the full width at half maximum (FWHM) decreases). At 
large energies it is greater for case A and C (almost co- 
incide) because of the larger distortion values. From this 
analysis it is revealed that the tunnel conductance can 
be improved in sp^-rich carbon for model C and also in 
sp'^-rich carbon for model B. Whereas for models A the 
conductance can vary as intermediately between model 
B and C. Although the tunnel resistance increases ini- 
tially in the low energy region it can be controlled for a 
wide band gap carbon structure. In practice the size of 
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FIG. 3. T{E) vs. E for the following sp^ bond concentration: 
(a) 72 % ,(b) 42 % ,(c) 24 %, and (d) 7.8 %. In subfigures 
(c) and (d) 3 curves (in color online) correspond to 3 cases 
of hopping disorder parameter 5 dependance on E-^ as shown 
in Fig. [IJb). The tight-binding parameters are the same as 
that in Fig. |l(c) The peaks for Model A showed a broader 
feature than that of case B and C, shown in the insets of Fig. 
(c) and (d). 
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FIG. 4. (a) Localization length (lower curves) and resistance 
(upper curves) calculated at Eres peaks, which are determined 
by the factor of disorder as well as the energy distance from 
the sp"^ structure band edge. The inset shows Lioc at low 
energies, (b) Variation of FWHM for the transmission peaks 
-Eres. 3 curves (in color online) correspond to cases A, B 
and C of hopping disorder parameter S dependance on Et^ as 
shown in Fig[2Ib). 



the sp^ clusters can be controlled to some extent through 
processes such as anneahng or irradiation hence the band 
gap can be controlled to some extent. Nitrogen doping 
(discussed belovif) is also a possibility. These results ex- 
plain the initial increase of resistance of tetrahedral a — C 
{ta — C) films incorporated with a small amount of nitro- 
gen followed by a decrease for high nitrogen (or sp^ — C) 
concentration [isj . 

A comparative study of the FWHM for cases A, B and 
C showed a significant difference not only in their abso- 
lute values but also in the their trends as a function of 
energy. For case A, the FWHM increases from low en- 
ergy followed by a saturation at energies above 2.5 eV. 
For case B, the FWHM decreases beyond 1.5 eV (due to 
decrease of disorder). For case C, the FWHM increases 
continuously above 2.5 eV. The other factors which effect 
the FWHM are the incoming energy Eres and the width 
of the sp^ "barrier" . A higher energy and smaller width 
increase the peak width due to the higher transparency 
of the structure for an incoming wave. The value of the 
FWHM, characterizing an average quasi-bound state life- 
time at a particular energy, approaches saturation in the 
high sp'^ limit. The width of the resonant peaks is propor- 
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tional to the characteristic time of resonant states there- 
fore as the disorder increases, the FWHM of resonant 
peaks and consequently the characteristic time of reso- 
nant states also increases. At the same time the ampli- 
tude of transmission maxima decreases rapidly. The con- 
ductance crossover occurs as with further decrease of the 
sp^ bond concentration the resonant peaks shift towards 
the mobility edge and the characteristic time therefore 
decreases with the resistance decreasing accordingly. The 
FWHM of the resonant peaks saturate at very low sp^ 
bond concentration but the resistance does not saturate 
as the amplitude of the resonance peaks increases as the 
peaks approach the mobility edge. These results clearly 
establish the specific effect of the non-linear model A, 
which shows the tuneable speed of carbon devices from 
the slow to fast regime as the sp^ % increase in these 
structures. 

The structural change for case A is also reflected in 
the LDOS accompanied by the peak position change and 
broadening. Fig. ([5]) shows the LDOS at the central cite 
for structures with different sp^ bond concentrations. At 
a high concentration (Fig. [5] (a)) there are several high 
LDOS peaks inside E^. As the concentration decreases, 
LDOS peaks decrease being comparable in amplitude to 
LDOS corresponding to the sp'^ bands (Fig. [5](b)). Fi- 
nally, when the sp^ bond concentration is less than 50 % 
the LDOS peaks inside become negligible compared 
to LDOS outside (Fig. [S] (c) and (d)). Also, the LDOS 
peaks move from the zero energy position with the con- 
jugated phase decrease, effectively opening the gap. The 
broad features of LDOS suggest the possibility for hop- 
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FIG. 5. LDOS (arbitrary units) for model A at the structure 
central cites for (a) 72%, (b) 62%, (c) 42% sp^ , and (d) 24% 
sp^ bonds (solid red curve) compared to that of pure sp^ 
structure (dotted black curve). 



FIG. 6. T{E) vs. E for the (a) 72 % ,(b) 42 % ,(c) 24 %, and 
(d) 7.8 % bond for the topological disorder considering model 
A. The energy axis is drawn relative to the Fermi energy. 



ping transport in these structures, which is commonly 
observed in diamond-like carbon films [3, Ull, [l^l • On the 
other hand, the high LDOS filling the sp^-gap region ex- 
plains the metallic conduction in graphitic carbon films. 
Hence the observed conductance crossover in a — C films 
with the increase of sp^C % (or N %) can be explained 
by the change of the characteristic time of the resonant 
states [25|. 

Besides the structural disorder we attempt to find the 
role of topological disorder on the transport properties 
of the structure. We know that topological disorder is 
more prominent in two-dimensional structures e.g. large 
sp^ clusters however, the effects can also be simulated in 
quasi-one dimensional (small sp^ clusters) structures. In 
the topologically disordered case for carbon superlattices 
or heterostructures, non-uniform sp^ — C cluster sizes 
were distributed throughout the system in contrast to 
the non-topologically disordered case where all sp^ clus- 
ters have the same size. Thus, it corresponds to so-called 
Anderson vertical disorder. The disorder parameter was 
determined using model A. As it can be seen from the 
Fig. [SI the topological disorder results in the decrease of 
the peak height and increase of the peak width and has a 
greater effect on the peaks located further from the Fermi 
level. The effect is more significant for structures with 
lower sp^ concentration, so the relative change in the qua- 
sibound energy levels is high. In addition to that, new 
minor peaks appear in when the difference in energy lev- 
els induced by topological disorder is larger than the hop- 
ping structural disorder parameter. Topological disorder 
modifies the resonant tunneling conditions, changing the 
resonant energy conditions of consecutive wells result- 
ing in a mismatch of the resonant levels. The amplitude 
of transmission peaks at resonant energies therefore de- 
creases and the resonant peaks spread out with energy as 
resonance no longer occurs at a sharply defined energy. 
A resonance peak splits in subpeaks corresponding to the 
cluster number (4 in this work), where topological charac- 
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teristic energy splitting Setop is proportional to ^ , where 
w and 6w are the average cluster size and its characteris- 
tic variation. When Setop is larger than the characteristic 
structural disorder variation determined via deformation 
potential as Ssstr — jSl, where 7 is the deformation po- 
tential and SI is the bond length variation, the subpeaks 
are distinguishable. With increasing structural disorder 
the subpeaks start to spread out and finally overlap when 
Sstop becomes smaller than Ssgtr with the elimination of 
the distinct subpeaks features due to the fact that the 
number of clusters is much smaller than the number of 
sites. The resonant peaks can therefore be tuned by dis- 
order of both types which can be controlled to some ex- 
tent experimentally through annealing or irradiation of 
the films fsd]. Based on the bond length distribution in- 
cluded in the model, different ring structures e.g. five fold 
and seven fold symmetry, additional so-called topological 
disorder can be included subject to further study. 

Conclusion: In this Letter, the importance of the non- 
linear disorder vs. energy model (case A), an interme- 
diate to uniform increase and decrease of disorder, is es- 
tablished. The theoretical analysis of this experimentally 
supported model discovers the possibility for a relative in- 
crease of the tunnel conductance in sp^C-rich carbon due 



to the appearance of resonant transmission peaks close to 
the bands. Nitrogen doping, for example, could be used 
to vary the micro-structure of a — C : films as nitrogen 
incorporation increases the sp^ — C ratio fis', 'so'] . The 
nitrogen concentration also infiuences the disorder [3l|. 
Nitrogen doping is therefore not the same as conventional 
doping however it can still be exploited to modify the 
conductivity of a — C films [18]. In this regard the result 
of the present study can be extended for doped carbon 
structures. A detailed study of nitrogen incorporation in 
a — C structures has recentlybeen carried out based on a 
model related to this work [32| . The lifetime of electrons 
can be nearly constant over a large range of sp'^C % in 
carbon structures having a wide bandgap energy. Since 
the intensity and position of the resonant peaks can be 
tuned by the sp^C % to sp^C % ratio and associated 
disorder under this model, we can effectively predict the 
nano-electronic device properties in a wide range of un- 
doped and doped carbon structures. This model can also 
explain the origin of observed resonant features in multi- 
layered carbon systems f25j. 
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